Non-adiabaticity and improved back-reaction 

Curtis Asplund t and David Berenstein tt 
' Department of Physics, University of California at Santa Barbara, CA 93106 
+ Institute for Advanced Study, School of Natural Science, Princeton, NJ 08540 

We treat quantum back-reaction in time dependent processes for quantum field 
theory in various simplified models. The first example is a harmonic oscillator whose 
' frequency depends on a second quantum variable x. Beginning with a classical 

^ ■ analysis, we show how using a particular canonical transformation the system can 

ry^ be described by an improved adiabatic expansion with a velocity dependent force 

for x. We find an instability at a critical velocity that prevents integrating out the 
oscillator degree of freedom in the new variables. We extend this calculation to the 
quantum system and to field theory and describe how to study fermions with similar 
techniques. Finally, we set up a model with an abrupt change in the oscillator whose 
quantum mechanics can be solved exactly so that one can study the effects of back- 

> 

■ reaction of a fully non-adiabatic change in a controlled setting. We comment on 

■ applications of these general results to the physics of D-branes, inflation, and black 



o 
o 
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holes in AdS/CFT. 
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I. INTRODUCTION 

Many problems in quantum field theory that involve the time evolution of a quantum system are very 
hard to solve. For example, thermalization of a system given some initial condition is such a problem and 
is crucial for understanding experiments like RHIC. Current methods in field theory do not explain the 
fast thermalization seen in this experiment [l]. 

As an approach to the problem of thermalization, one can study the excitation of degrees of freedom due 
to the time dependence of some other variables. Such problems typically arise in the study of cosmology, 
in particular the problem of reheating [2]. Similar problems arise in string theory when one studies the 
collisions of D-branes and string creation between D-branes. 

These problems usually start with some degrees of freedom that are in motion and some other degrees 
of freedom that become excited due to this motion. The first set (those that are initially in motion) are 
usually called the moduli fields. Let us call them 0. During this motion there are other degrees of freedom, 
e.g. the modes of a massive field, that are sensitive to the motion of the moduli. Let us call these the 
heavy degrees of freedom. The heavy degrees of freedom have large associated frequencies w(0), so these 
degrees of freedom are fast. In most field theory setups these can be thought of as harmonic oscillators. 

Let us suppose they all start in the ground state. If the moduli motion is classical and the velocity 
small, then the adiabatic theorem guarantees that they stay in the ground state |3|]. However, if the 
velocity of the moduli fields is large enough or if oj(4>) small enough, there can be transitions to excited 
states. This results in particle production in the field theory case. This describes the effect of particle 
production in cosmological evolution (see [4] for a pedagogical example). In most cases, the moduli field 
motions are treated as a background on which one does computations for the heavy degrees of freedom, 
which simplify to a set of decoupled harmonic oscillators with ui((j)(t)) as time dependent frequencies. 

In such situations, the full calculation in the background field consists of a Bogoliubov transformation 
between the Fock space of heavy states in the initial time and the Fock space of heavy states in the final 
time [5]. Here the particle production is considered to be small, so that interactions between the particles 
that are produced and the effects of these particles on the background moduli can be ignored. With some 
extra effort, dissipation can be added to the background moduli to account for some of the energy lost to 
the particles that are produced. 

We are interested in a more systematic treatment of the production of particles with the back-reaction 
of the moduli included. We want to improve on the adiabatic approximation and to understand how it 
breaks down. This is particularly important in setups where one expects strong back-reaction due to the 
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excitation of the heavy degrees of freedom, say due to fast thermalization. This mechanism is believed to 
operate in the problem of black hole formation in the dual CFT of AdS / CFT setups [6] and this is one 

of the main motivations for this work. The black hole formation process has been well studied on the 

I I 

gravity side [7[ (and references therein), but the problem of the time evolution and thermalization of the 
quantum field theory is largely unsolved. However, see js] for a solvable toy model that addresses some of 
these issues. 

For all of these problems we want to study what happens once we take into account this particle 
production mechanism in the dynamics and include it in the evolution of the moduli fields. It is important 
to approach the problem with caution, for if we think of particle production as a measurement, the different 
particle production outcomes could decohere and each such state would then evolve independently of the 
others. With this in mind we do not average over the heavy degrees of freedom, as this might lead to 
wrong results. 

In this paper we begin to systematically explore these effects in simple toy model systems. We will 
have one modulus (light) field and one heavy field whose mass depends on the modulus field, treating 
the cases of heavy bosons and fermions separately. We will find it more convenient to work in the 
Hamiltonian formalism. The process of integrating out degrees of freedom in a background field is more 
common in the Lagrangian formalism in field theory applications. This same procedure of integrating out 
fast degrees of freedom is performed in the Born-Oppenheimer approximation in the study of molecules 
or solids, where there is also extensive work on higher-order corrections [91]. A lot of our work is to 
rewrite the Hamiltonian in variables that make this procedure more transparent and that permit one to 
study systematic improvements to the adiabatic approximation in a manner and setup appropriate to the 
applications we have in mind. 

In the next section we will write a Hamiltonian for the light and heavy fields in a basis that diagonalizes 
the heavy field degrees of freedom for each value of the modulus, considering first the classical then the 
quantum problem. Here the appropriate adiabatic approximation is the Born-Oppenheimer approximation 
where the heavy fields are integrated out. We keep a parameter that controls how much back- reaction is 
present, the classical mass of the modulus degree of freedom. At high mass and fixed velocity, the kinetic 
term of this degree of freedom dominates the energy, so the force due to the heavy fields not being in their 
ground state has a negligible effect. That is, the degree of freedom has a lot of inertia. For low mass, 
any small modification of the potential has a large effect, because the kinetic energy is small. In the third 
section we discuss what happens in the case of heavy fermions. In the fourth and final section we set up 
a simplified model with an abrupt change in the oscillator degree of freedom whose quantum mechanics 
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can be solved exactly and we investigate and interpret the solutions. 



II. THE MINIMAL MODEL 

We begin with a simple quantum system with two coupled degrees of freedom, one of which can 
interpreted as a rolling modulus field. This is a minimal model for studying back-reaction in the regime 
beyond which the adiabatic approximation breaks down: we need one adiabatic degree of freedom that 
turns non-adiabatic and we also need another degree of freedom on which to back- react. Our model is 
specified by the following quantum Hamiltonian 

This is the Hamiltonian for a particle in two dimensions with potential U (x, y) = \Vt 2 (x)y 2 (in units where 
h = 1). We normalize the mass of y to unity but keep the mass parameter m x of x. The function Q controls 
the curvature of the quadratic potential well in the y direction. A similar toy model was considered in 
the problem of D-brane scattering [l0|, where Q(x) oc \x\. A classical analysis of this simple model can be 



found in 
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The case of a harmonic oscillator in the ^/-direction with angular frequency u is just Q(x) = to. In that 
case the Schrodinger equation may be easily solved (via separation of variables), yielding the spectrum of 
eigenfunctions with wavefunctions 

j> k>n (x, y) = e kx N n R n {uj x l 2 y)e-^l 2 (2) 



where N n = (w/vr) 1 / 4 (l/Vn\¥) is a normalization constant and H n is the n th Hermite polynomial. Here 
k G K. and n G N. The corresponding energies are 

Ek, n = k 2 /2m x + uj(n + 1/2). (3) 

In a classical setup we can always set y = p y = for any Q(x) and get a solution of the equations 
of motion. The x variable then behaves as a free particle and the solution to the equations of motion is 
x = vt, where v is some velocity. Thus x essentially measures the time. We can even consider a classical 
potential for x if we want to, but the point of our paper is to understand y and how this degree of freedom 
back-reacts onto the motion in x, especially in the quantum system, so we will not study the dependence 
on a potential in x. 
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Assuming x = vt, we can then study infinitesimal variations around this solution with y and p y small. As 
such, the y and p y infinitesimal motions represent a harmonic oscillator with a time dependent frequency, 
given by u(t) = Q(vt). The motion in these variables is considered adiabatic if 

w n' , ... 

This can be achieved with small velocity, large Q or a small gradient of Q. But if we increase v enough, 
we can always go to a velocity regime where the adiabatic approximation breaks down. We will see that 
it is exactly in this regime that we can not ignore quantum corrections. 

We would like to understand the breakdown of adiabaticity of this system, in the form of a variation 
of Q(x). That is, we would ultimately like to have an expansion (of, say, the Hamiltonian) in terms of 
derivatives of Q that lets us systematically correct adiabatic or Born-Oppenheimer approximations and 
lets us know when such approximations break down entirely. Remember also that in a quantum system 
Q(x) gets replaced by an operator, so the use of x as time becomes more of an issue. 

In standard semiclassical treatments, one would treat x classically but solve for y quantum mechanically. 
Motion in x can lead to particle production (excitation) in the y direction. The equations of motion of x 
can then be corrected by averaging over the particle production of y. However, as we mentioned in the 
introduction, this averaging procedure can be problematic for solving the real quantum problem accurately, 
so we do not do this here. 



We can vividly see the breakdown of adiabaticity by first considering the system classically. The classical 
limit of this quantum system is given by 



The appearance of y in the solutions always as uj x l 2 y suggests that a change of variables y = VL^^y may 
be useful in understanding the essential degrees of freedom in the system. So let us consider the system 
in terms of the variables y and x = x. We have that 



A. Breakdown of adiabaticity 




(5) 




(6) 



(7) 
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These can then be substituted into the Hamiltonian. There is also a change of measure dxdy = Q^-^dxdy, 
and one has to be careful about this. We will address this in the quantum case in the next subsection. 
To implement this change of variables in phase space, we make the canonical transformation 

x := x (8) 

y:=n l l\x)y (9) 
1 Q'(x) , , 

n '= Px ~ 2W) yPv (10) 

Canonicity of this transformation can be checked from the invariance of the Poisson brackets. Notice that 
the change of variables in y,p y — > y,Py is obviously a change of scale that is x dependent and that their 
Poisson brackets are retained. Also, the Poisson bracket with x = x vanishes as before. However, p x does 
not Poisson-commute with y,Py and it gets corrected. This can be guessed by noticing that the bracket of 
p x with y is proportional to y, and the bracket of yp y with y is also proportional to y. Therefore we can 
try to cancel these against each other, getting the result above. There will be a similar situation when we 
deal with fermions. 

The Hamiltonian in these new variables is 

2 „.,„...,. a . n V 

2m r 2 y 2m r \l 8m r \ \l 



where Q = Q(x) and Q' = dQ(x)/dx. Notice that the Hamiltonian receives both quadratic and non- 
quadratic corrections in the y,Py variables. The coefficient of the former is given by 

Ps n* in 

(13) 



2m x Q 2fi' 

the latter expression valid in the y = p y = 0, Pxjvn x = v regime. 

From this Hamiltonian we can already see that there will be a qualitative change in the behavior of 
the y degree of freedom as Q' differs from zero. For constant Q, the Hamiltonian for y and Py is simply a 
classical harmonic oscillator, with phase portrait as shown in Figure [TJ 

As vfl'/fl becomes non-trivial, the third and fourth terms in ( Tl2|) become important, changing the 
orbits from circular first to elliptical and then to hyperbolic. At that point the origin becomes an unstable 
fixed point as shown in Figure [2j It is important to remember that here we are taking slices of phase space 
in the x, y variables and not the original x, y variables of the system. 
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FIG. 1. The phase portrait of y and py for constant SI 




FIG. 2. The phase portrait of y and p y for non-negligible f2'/f2 

We can be more precise. The behavior of the Hamiltonian vector field near the origin of phase space 
of y is determined by the second degree terms in y and p y in ( f!2l) . The stability of the solutions is set by 
the signs of the eigenvalues of the quadratic form 



ft 



(14) 



Q/2 (tt'/tt) Pi /4:m a 

(the overall factor of 2 is added for later convenience). The condition that the origin is stable is that the 
determinant of this matrix is positive (Q is positive in our conventions). This can be seen by diagonalizing, 
which puts the Hamiltonian in the form of a simple harmonic oscillator and the sign of the determinant 
then corresponds to the sign of the quadratic potential. The condition for stability then gives us 

fi' p- x 



W2m x <l - (15) 
This is analogous to the non-adiabaticity condition uo/uo 2 < 1, often used as a threshold for non-adiabatic 
behavior or particle production. Here it sets the threshold for treating y,Py as a harmonic oscillator versus 
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an oscillator with an inverted potential (this requires rotations of the variables to make it explicit). That 
is, when we have a classical motion in x the y degree of freedom becomes unstable exactly when 

(W ) 

and we expect to have corrections of order one in the y degree of freedom, when u/u 2 ~ 1. Indeed, we 
can improve the adiabatic approximation with the diagonalized quadratic action in y,Py. Recalling that 
the canonical form of the harmonic oscillator is the diagonal form p 2 /2m + mu 2 x 2 /2, the diagonalization 
gives us an oscillator with a (position and momentum dependent) frequency 



^Vd^/ff-jLg) 2 (it) 

If we consider a semiclassical treatment with x = vt to quadratic order in y, the zero point energy of 
such an oscillator is then 



1 1 
-co = - 

2 2 



\ 




+ . . . 



where the last expression holds for small Q/Q 2 and we have used Px/m x ~ v and Q = Q'v. We see that 
in this way it is possible to generate velocity dependence in the low-energy effective Hamlitonian for the 
moduli degrees of freedom, in this case \{p\/m x + Cj). Notice that there is apparently already an x- velocity 
dependence in the term in (fl2|) linear in p x . However, when integrating out y and Py one finds that to first 
order in Q this term has no contribution since ypy averages to zero over the motion in this regime. Thus 
this term should contribute only at second or higher order, consistent with our result above. 

Velocity dependent terms like this have been argued to give rise to gravitational interactions in matrix 
theory as a dual description of M-theory in 11-D 12j. Also notice that the square root formula is remi- 
niscent of a Dirac-Born-Infeld Hamiltonian for D-branes, but the sign of the velocity term in the square 
root is wrong for that comparison. As we mention in the conclusion, a full investigation of these issues 
and connections is reserved for future work. 

When the origin in the y,p y plane becomes unstable, the fact that in a quantum system a state 
occupies some finite area of phase space implies that the quantum state will spread along the trajectory 
of instability. At that stage we cannot integrate out the y variable any longer and it has to participate in 
the full dynamics. This is a nice semiclassical way to explain that there must be particle production at 
large v compared to the case of v near zero. In the following sections we will turn to the quantum system 
and explain how to treat the adiabatic approximation and its breakdown more carefully. 
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We can interpret the effect on the x degree of freedom as a momentum dependent force 0. Indeed, let 
us generalize the above result and consider a case where the degree of freedom x is some rolling modulus 
field and we have field theory degrees of freedom r/k labeled by their momenta in a <i-dimensional box of 
volume V. The corresponding field theory would live in d + 1 dimensions. 

The effect we consider will just sum over the modes. For each of them we will find the usual relativistic 
dispersion relation, where Q(x) plays the role of the mass. The frequencies for the mode labeled by the 
momentum k is given by 



tt k = + k 2 (19) 

and thus 

* = vPTP' (20) 



We find this way that the modified frequencies ( fTTj) are given by 



^k 



^ + k 2 -\[-^-\. (21) 



The zero point energy contribution (after subtraction of the term at zero velocity in order to get a finite 
answer) gives us an additional contribution 



i£ = * = = £5 .«S-i | - 1 -n^-^n A nn r + .... (22) 

k k \ \ / k 

The sum is convergent if the number of spatial dimensions is less than five, d < 5. The effective term in 
the hamiltonian, after integrating over all modes, is proportional to 

where V is the volume of the d dimensional box: it fixes the dimensional analysis in the equation. We 
find this way a finite effect per unit volume, with a fixed sign. In three spatial dimensions d = 3 we notice 
that the correction only depends on the gradient of the mass and the velocity of the modulus field. This 
can be generalized to multi-field moduli, giving essentially the same result, where now Q = VjfivTj/m and 
7Tj is the canonical conjugate of x l . This first term in the correction can be interpreted then as a modified 
metric on the moduli space (a sigma model would be written with a Hamiltonian given by H = g % ^ r n i i\j\ 
It has a definite sign, making distances longer along the gradient of Q (the sign is negative above because 
it is associated with the inverse metric) and giving rise to a curved metric on the moduli space. 
1 Such terms are usually computed as a correction to the action in the path integral formulation around a classical solution 
of the x motion. This is exemplified in D-brane setups in 13| , where one also sees the appearance of particle production 



from the imaginary part of the phase shift. However, no back-reaction is taken into account. 
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B. Quantum adiabatic approximation and corrections 

Now that we have described the classical theory with a brief excursion into a simplified description 
of the quantization, we turn our attention to dealing with the fully quantum problem. Within the adia- 
batic regime where (1151) is satisfied we are motivated by the above canonical transformation to consider 
wavefunctions of the form 

IK*. V) = E /»(*) (—) V4 ^=H n {^\x)y)e-y^\ (24) 

n=0 \ / v 

and so specify a state by the set of functions {f n (x)}. This is a change of basis or representation in the 
Hilbert space of states. There is no loss of generality from writing the wave function in this fashion and 
the following results are still exact. 

Wavefunctions of this form have the nice property of being energy eigenstates with respect to the y parts 
of the Hamiltonian (at each x). This setup and the following analysis is exactly the Born-Oppenheimer 
method of solving the dynamics of quantum systems. One first divides the variables into slow and fast 
degrees of freedom. One then solves the exact dynamics of the fast degrees of freedom when the slow 
degrees of freedom are frozen and compute the Hamiltonian in the new basis. The further assumption to 
simplify the problem (the Born-Oppenheimer approximation) is that no transitions occur between different 
levels as the slow degrees of freedom move, mainly that the motion is adiabatic. This approximation has 



been studied extensively. We would like to draw particular attention to the works of M . Berry [14| where 
this setup was explored systematically. Apart from the Berry phase connection, this generally leads to 
effective terms in the theory that can be attributed to effective electric fields and magnetic fields for the 
slow variables (see 15| and references therein) 3- Our results do not assume the adiabatic approximation, 
but we do perform the appropriate adiabatic change of variables. 

We have chosen the wave functions for y to be normalized to one at fixed x. This means that the change 
of measure from the x, y variables to the x = x, y = yQ 1 / 2 has been absorbed in the normalization of these 
states, which is why we encounter Vt(x) l l A sitting in front of the y wave functions. This implements the 
quantum change of variables so that the y part of the wave function has the correct behavior with respect 
to y variables in the number occupation basis for the oscillator of frequency Q(x). 

We now compute the matrix elements of the Hamiltonian in this basis of states. Non-adiabatic changes 
in Q will manifest as off-diagonal contributions to the Hamiltonian in this basis, which is what we're 
interested in. In a general energy eigenstate basis {jV'i)} we can ac t on a state \if)) = J^n^l^n) with 



Many of the standard results of D-brane interactions 



due to integrating open strings can be rephrased in this language. 
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H to obtain H\ijj) = E n c n \i/) n ). That is, we obtain a new state with coefficients weighted by the 
corresponding energies. Here we do an x dependent version of this, which will tell us how the Hamiltonian 
mixes states of the form (|24"1) . 

We write the wavefunction of H\ip) in the same form as ( 124|) . with a new set of functions {g m (x)}. 
Then we can extract a matrix representation of the Hamiltonian by using orthonormality and integrating 
out the y variable: 

9m{x) = y^ j H mn (x)f n (x) (25) 

n 

where (no summation on n) 



H m nfn 



n 



X 



y/m\2 m Vn!2" \vr 
Id 2 Id 2 



1/4 



2m T dx 2 



2 dy 2 2 V 



J H m (fl 1/2 (x)y) e-y 2 ^ 2 
fi(a;)^ 1/4 



(26) 



fn(x) 



IT 



H n (Q 1 / 2 (x)y) e 



dy. 



One term in the diagonal part of this matrix can be computed immediately, since for fixed x the factors 
H n and the exponential are an eigenfunction of —d 2 /2 + Q 2 y 2 /2, with eigenvalue Q(n + 1/2). So along 
with the orthogonality of the Hermite polynomials we have 

Hmnfn = fm^(m + l/2)5 mn + (integral with -dl/2m x ). (27) 

The computation of the remaining integral is straightforward but long, as the second derivative gives 
rise to 16 terms. But each can be evaluated using the orthogonality and recursive properties of the Hermite 
polynomials. The result is 



Hry 



Jn = ~ ^ — V(™ +l)(m + 2)(m + 3)(m + 4) /, 
62m r 



12 



y/(m + l)(m + 2) 



5m, 



+ 



+ 



f m Q(m + 1/2) 
^Jm(m — 1) 



I [^-^\+2f ^ 



C Q' 2 m 2 + m+l 

i" 1 1 



(I 



m(n—2) 



2m T 



Q? 16m T 



8„ 



8m x 



n" n ,2 \ 



fm-2 -TZ T^T + 2 /m-2 7T 



n n 2 J 



5, 



m(n+2) 



ZQ^VMjn ~ X )( m - 2 )( m - 3 ) /w-4 7^-^m(n+4)- 



(28) 



Here we have an explicit expansion of H in terms of the derivatives of f2 □. The adiabatic approximation 
in this context is to ignore all terms in the Hamiltonian that depend on derivatives of fl Notice that these 
3 This is the same expansion as in the Born-Oppenheimer method, where the fast degrees of freedom are oscillators rather 
than electrons in a molecule or material. The off-diagonal terms correspond to the "vibronic couplings" or "non-adiabatic 
coupling terms" in that context 9j. Our results in this section are a specialization of the Born-Oppenheimer method to a 
case where the fast degrees of freedom are oscillators, which can be interpreted as modes of a massive quantum field. 
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are the only terms that can change the level of the oscillator, so in the adiabatic approximation we end 
up projecting onto the components of H that lie on the diagonal H mn oc 5 mn , giving the result 



The system splits into an infinite number of one dimensional problems for a particle on a line with potentials 



These define what are called the potential energy surfaces, with a conical intersection at x = 0, in the 
Born-Oppenheimer literature. We see here the usual energy of a harmonic oscillator in the y variable 
including the zero point energy. 

We can now expand on the terms that are ignored in an adiabatic approximation as a perturbation of 
the Hamiltonian. The terms that involve x-derivatives of the wave function coefficients /„ only change 
the level by ±2 units (in particular, they keep the even/odd splitting of the y wave functions). The terms 
that do not involve x-derivatives of the /„ involve two x-derivatives acting on polynomials of Q. Only the 
first set of terms, involving an f' w are momentum dependent. If we compare these with our classical result 
given in equation (I12p . we find similarities in the description. In terms of lowering and raising operators, 
we usually have that y oc a + a) , while pg oc a — a). The term ypy is then of schematic form a 2 — (e^) 2 , up to 
normal ordering ambiguities. A careful analysis shows that the off-diagonal terms that involve derivatives 
of the form f' m exactly match the harmonic oscillator algebra coefficients. Thus, we can state that in the 
quantum theory, the momentum dependent corrections depend on the operator 



For semiclassical x motion, p x /m x ~ v and if we take into account only these corrections, we can obtain 
a power series in v for the ground state. This arises from the self-consistent assumption /„ ~ 0(v n ), 
beginning with /o ~ 0(1) (but this can be generalized to other states). The corrections to the energy 
etc, will then always be even powers of the velocity. These terms survive in the m x — > oo, v fixed limit. 
This limit is the classical limit for x. At large mass for the motion in x, fixed f2(x), there is essentially 
no back-reaction. This is because the system is dominated by the kinetic energy E ~ \m x v 2 and Q is a 
parametrically small perturbation. Also notice that the terms that do not involve momentum are those 
that correspond to the square ijjpy) 2 and these are suppressed in the large m x limit. 

We can also work in a level truncation approximation (again removing the terms with no /'), where we 




(29) 




(30) 




(31) 
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keep only the modes n = 0, n = 2. The Hamiltonian in two components will then be of the form 

■(;:)= (t .;) «s (i, :))(;:) « 

If we set /o — Coexp(ikx) and /2 — C2exp(ifcx) to check how important the off-diagonal terms are (again 
in a limit where x is essentially classical), we have a level splitting of 2Q compared with the off-diagonal 
component that is y/2Q! /Vl-^-. 

Diagonalizing the above Hamiltonian we find that the lowest energy level for the two-level system is 

We see that to second order in the Taylor series in Cl we find a match with our semiclassical analysis in 
(TTSJ), making the naive semiclassical result more plausible. 

One should notice that the sign of the term that includes Q 2 is negative. This is expected because 
the perturbation is off-diagonal and under this condition in second order perturbation theory in quantum 
mechanics the correction to the energy of the ground state is non-positive. 

If we go beyond the adiabatic regime and pass through a region where the y variable is unstable, the 
semiclassical result ( II 8p becomes purely imaginary in this region. The integrated effect of this term can 
be interpreted as an imaginary contribution to the action for a persistence amplitude in the ground state, 
indicating that the system ends up in a state different than the ground state with some finite probability. 

It's useful to estimate the probability that we go from the level with zero occupation number to the 
level with occupation number two in the regime where there is no back-reaction on x. Remember that 
having no back-reaction on x is equivalent to the limit where we send the mass of the x degrees of freedom 
to infinity, keeping the velocity fixed. Also, x(t) becomes a classical trajectory: at finite velocity the 
momentum conjugate to x is very large and we can localize the wavepacket in the x direction to arbitrary 
precision without having to take into account the finite time broadening of the wave packet. 

In such a system, we want to understand the dynamics of the y degree of freedom. For simplicity, we can 
truncate the system to two levels as above, in the regime where the off-diagonal terms are parametrically 
small (the truncation to these levels is only valid in this regime). Then Q(x) ~ Q(x(t)) and the system is 
a two level system with effective Hamiltonian given by 

(2VL + \VL \/2i%\ 
H{t)=\ * \ 4 • (34) 
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In such a system the off-diagonal terms constitute the perturbation Hamiltonian. The difference in energy 
between the excited state and the unexcited state means the system will be subject to Rabi Oscillations. 
We will see that we get a small average occupation number if the off-diagonal term is parametrically small, 
i.e. the system is near adiabatic. 

The amplitude is given to first order in time dependent perturbation theory by 



A)->2 = 7^7= J dt ^ exp ds 2ft(ar(s))J . 



(35) 



To estimate, we take Q large, Q small and both approximately constant over a time T. This is then 
roughly given by 

A^ 2 ^-^~J rftexp(2z^) = ^=^exp(^T)sin(fiT). (36) 
The probability of finding the state in the excited configuration will then be given by 

P(T)^sin 2 (fiT), (37) 

where rj = Cl/Q 2 is the adiabaticity parameter. 

This averages over time to r] 2 /lQ. So we find that the probability of being excited grows like the velocity 
squared, but stays small in the adiabatic regime r] <C 1. Notice that this is a short time estimate (we 
assumed Q approximately constant). We can do a similar calculation if we already begin in an excited state 
at level m and consider transitions m — > m ± 2. We find that the amplitude to excitation is slightly larger 
than to de-excitation because of the harmonic oscillator algebra. This is similar to stimulated emission of 
radiation: we get signal amplification. This is typical for this class of problems [5(. 

Notice that this is very different from particle production for an infinitely long time process. For many 
of these the amplitude over an infinite regime becomes non-perturbatively suppressed in the velocity for 
small velocity, so long as there is a minimum frequency in Q(t). This is typical of D-brane scattering with 



finite impact parameter. This is nicely explained in 



221. 



The main reason for extra suppressions at large times is the following. At constant velocity for x we 
have 

tt(x(t)) = Q(vt). (38) 
Then Cl/Q 2 = vQ' /Q 2 and the phase in the exponential is given by 

t 1 px=vt 

ds 2Q(vs) = - du 2Q(u). (39) 

v Jo 
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The transition amplitude takes the form 

= 571 L ix ~^ exp [v I iu 2Sl(u) ) • (40) 

which oscillates very quickly when we send v — > and the behavior of the integral is similar to that of 
the Fourier transform of 0! /Vt 2 at high frequency. We have in mind changing variables to 7 = Vt(u)du. 
The resulting integral receives contributions from the poles in the complex plane where ^(#(7)) = and 
this can be very different from integrals at finite time where we cannot deform the contour. 

Notice that this number is generally much smaller than A ^ 2 for finite time, where the result is polyno- 
mial in the velocity. The result that is polynomial in the velocity contributes to back-reaction as a power 
series in v, matching the type of expansion we saw before and what is expected from direct calculations 
of velocity dependent forces. 



III. FERMIONS 



We now discuss fermions for completeness, both at the semiclassical level and at the quantum level. We 
consider systems where there is a fermion parity operator (— 1) F , with F counting the number of fermions, 
that describes whether a state is fermionic or bosonic and is preserved by the action of the Hamiltonian. 
Since we wish to study a system that can mix with a bosonic system, we need at least two different 
fermions, giving us at least four states: two of even parity, two of odd parity. The even states can mix 
with each other, but not with those of odd parity. To follow the discussion of bosons above we will need 
to use the Hamiltonian formalism for fermionic (anticommuting) variables, as described below. We begin 
with a single fermion oscillator then consider the possibilities when additional fermions are present. 

A fermion oscillator is defined with a pair of Grassmann variables 9 1 ' 2 with Lagrangian given by 

L = X - [i9 x 9 x + i9 2 9 2 - Mij&Qi} . (41) 

The matrix My is antisymmetric, and the classical variables anticommute with each other: {9 l ,9^} = 0. 
The 9 variables are real: 9 l * = 9 l . The factors of % in the kinetic terms are there to make the Lagrangian 
real (using the usual property of complex conjugation that also reverses the order of the variables). The 
Hamiltonian is given by ^M i j9' l 9 j and it is real if M is antisymmetric and hermitean. The (left) canonical 
conjugates of 9 1 ,9 2 are 717 = —i9 1 /2, tc 2 = —i9 2 /2, and the Poisson brackets are 



(42) 
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Here we are following the "superclassical" formalism and conventions of 171 ] . The equations of motion 
that follow from the Lagragian are 

V 

e 2 





(43) 



corresponding to an oscillator with frequency to 
the fermion oscillator see 18 ]. 



|M 12 | = V det M. For a detailed solution and more on 



Now we consider coupling many such fermions to a classical variable x, so that My(x) is an 2N x 2N 
antisymmetric, hermitean matrix that depends on x. The trick to solving the fermions is to rotate all the 
variables 9 l = R}^ into each other with an SO(2iV) transformation so that M takes a standard form 

/ M 12 {x) ...\ 

-M 12 (x) ... 

M{x) = R T MR = M u (x) ... ■ (44) 

-M 34 (x) 

V : • i •7 

The Poisson bracket of x with Q % vanishes, but the Poisson bracket of p x with 6 l does not, analogous to 
the classical bosonic system we first considered. Indeed, we find that 

{p x , 0%b = ~d x R)W = -d x R)(R7 l )i6 k . (45) 

Now, we can correct p x by defining 

p x = p x + SfaQi (46) 

so that the Poisson bracket with 9 l vanishes. We find that S = R~ 1 d x R. With x = x and 7fj = (R~ 1 ) J i 7rj 
we have a (super) canonical transformation. 

The quantization of the fermion oscillator is somewhat subtle, due to the fact that it is a constrained 
Hamiltonian system. In particular, it has the (second class) constraints iTi + {i/2)5ij6 : > = 0. These can be 
used to eliminate the momenta (on shell). Because of the constraints naive canonical quantization cannot 
be applied direeth, using the Po 1Ss o„ b r ac k et, but m ay be lm ple m euted via the Dhac bracket Q or the 
Peierls bracket [18(, leading to the (equal time) operator anticommutator 



{§\e j } = 8 ij . 



(47) 



From these we can define an annihilation operator a = -^(O 1 —i0 2 ) and its conjugate creation operator a\ 



in terms of which we can write the Hamiltonian: H 



ooia'a 



trt — i 



). The spectrum of H is {— to/2, +u/2} 
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and the non-degenerate Hilbert space is two-dimensional. If we have iV fermion oscillators then the Hilbert 
space is of dimension 2^. 

We can discuss corrections to the strictly adiabatic case that are made manifest by the above super- 
canonical transformation, in parallel with our discussion of the bosonic case. In the transformation x 
enters through the dependence of R(x). Remember that R is a matrix that block-diagonalizes M. If we 
have only one fermion, R is trivial and the state with the fermion occupied and the fermion unoccupied 
do not mix. This can also be seen from conservation of fermion parity symmetry. 

In general we get corrections in derivatives of R. For two fermions there are four states in the quantum 
theory, two with even parity and two with odd parity. The two with even parity form a two-level system 
with a similar Hamiltonian to the one we computed in (|32|) . except that the off-diagonal terms are not just 
(Q'/Q)v, but involve the components of the dR/dx matrix. These velocity corrections are qualitatively 
similar to those for bosons truncated to a two level system: after diagonalizing we will get a square root 
formula with momentum dependence for the zero point energy. Notably the correction to the zero point 
energy carries the same (negative) sign as in the bosonic case. 

If there is no mixing between different fermions because M is strictly block diagonal for all x, then we 
find that there is no fermionic contribution to the momentum- dependent forces, nor any fermion "particle 
production" unless M changes sign (this is what happens in the simplified model studied by IQj). Similar 
effects were found in 



IV. ABRUPT CHANGES 

We can also consider non-adiabatic situations. In general the non-adiabatic case is not solvable and one 
has to resort to simulations. Here we consider a simple example that is tractable so that we can extract 
lessons for more general cases. 

A particularly simple and extreme example involves an abruptly changing potential where 

f 1 : x < 

Q 2 (x) = { ~ (48) 

I u 2 : x > 0, 

with u> 2 > and w ^ 1. In this setup the wave functions on the left x < and right x > are solvable 
as they reduce to two separable degrees of freedom. The system reduces to solving a boundary condition. 
Notice that in this case we can always choose the mass of the degree of freedom x to be equal to one by 
rescaling. We expect both reflection and transmission from such a potential, which is depicted in figure O 
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FIG. 3. The abruptly changing potential 



The physical situation we are modeling is a particle coming in from x = — oo and being scattered by 
the change in the potential. Alternatively, this represents an object with two degrees of freedom, one of 
which represents some center of mass motion (the x variable) and another that represents some internal 
excitation (the y variable) whose structure changes at x — 0. 

To solve this problem, we solve the Schrodinger equation away from the origin, imposing boundary 
conditions on the solution appropriate to the physical situation described above. One is free to choose an 
incoming state. We impose that the incoming particle has momentum k and that it is in the ground state 
with respect to the y variable, hence it has energy E = k 2 /2 + 1/2. The classical limit is where the energy 
stored in the x degree of freedom is very large, expressed by a large k. 

Left of the origin (i.e. x < 0), we have the normalized solution 



We always take the positive branch of the square root, ensuring that for I < k 2 /2 we have a propagating 
mode to the left and for I > k 2 /2 we have a decaying mode to the left. If this were a field theory problem, 
we would need to match the energy with zero point subtraction, which would be the same since both of 
them are to the left. 




1=0 



(49) 



The sum over modes is the set of reflected waves with excitations of the y degree of freedom. 
Conservation of energy requires that E = k 2 /2 + 1/2 = kf/2 + I + 1/2 for each I, or 



h = Vk 2 - 21. 



(50) 
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In the positive x region we can similarly write the wave function as 

1 



iKs, y) = J2 C h Q -j=H h (V^/K^ W, (51) 



h=0 

,2 In , 1 /o _ 7V 2 



with conservation of energy here requiring .E = A; 2 /2 + 1/2 = k' h z /2 + u(li + 1/2), or 

fcjj = v / ^ 2 + l-w(2/ 1 + l), (52) 

again taking the positive branch so we have propagation or decay to the right. We choose no incoming 
waves from the right. 

In a field theory setup the zero point energy on the right and left would be different and in general 
it would be infinite. To match we would need to work out a renormalized model properly taking into 
account the finite parts of the zero point energy. Here we can model it with an offset in the conservation 
of energy if we want to. Only the conservation of energy equation for matching between left and right 
would require such an offset. In this case it would represent an additional contribution to the potential 
in the x direction given by a6{x). In what follows we will not consider this possibility further: the results 
are similar to those with a — 0. 

Our goal now is to calculate B[ and exactly. These coefficients represent an S-matrix type problem 
between many channels meeting at a boundary. We must impose continuity of and dip/dx at x — 0. To 
accomplish this, we need to express H n (y/ujy) in terms of H n (y) via a transformation matrix, U mn . This 
can be derived directly from the generating function for the Hermite polynomials. It can also be computed 
from the creation and annihilation algebras as well as the corresponding vacuum states in the two regions, 
related by a Bogoliubov transformation. We describe this latter calculation in Appendix 1, where we also 
describe some further properties of U mn . We find that 




1/4 '^V4 1 



71^ y/m\2 m Vn 



1 f°° 

= / H n (x)H m (^x)e- x2{1+ ^/ 2 dx (53) 

l\2 n J-oo 



2wV2 1 



l + w((m + ra)/2)! 

a k (-a) k+{m - n)/2 (3 n -' 2 \ (54) 



(m + n)/2 \ k( \k+{m-n)/2an-2k 



K k, k + (m — n)/2, n — 2k 
where the first factor in the sum is a trinomial coefficient ( 6 "J = an d 



1 — to , , 

a := 55 

1 + w 



Am 1 / 2 
1 + cj' 



(56) 
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In the formula we must also add that the integral is zero unless m and n have the same parity. By 
orthonormality of the Hermite polynomials (properly normalized and weighted), we may write 

1/4 



W\V4 1 



7T/ Vm!2 r 



/ 1 \ 1 1 



(y)e' 



(57) 



or in terms of the harmonic oscillator stationary wave functions 

'un 1 /* i 



u n (y,u) :-- 



H n (y/ux)e' 



-y 2 Lo/2 



we have 



(51 



(59) 



n=0 



Note the matrix U mn is a change of basis matrix between two real orthonormal bases, and as such is 
orthogonal. 

Now we can impose continuity of ip at x = 0. Approaching from the negative x region, we have 



V>(0, y) = Au (y, + B i u ^ V < 



(60) 



and from the positive x region 



^(0,y) = ^2C h u h (y,u) 

CO oo 

= ^2 C h^2u hl2 u l2 (y,l) 



l 1= l 2 =0 



C h U hh u h (y,l). 

h,h=0 



(61) 



Now we can match the coefficients of the ui in (160]) and (1611) . which we can write as a matrix equation: 



( A . . . ) + ( B B x . . . ) = ( C d . . . 



( U 00 U i . . . ^ 



(62) 



or more succinctly A + B = CU, with obvious matrix notation. 

Similarly, we can impose the continuity condition on the derivative dip/dx. Coming from the negative 
x region we have 



dip 



— (0,y) = ikAuoiy, 1) -^ihBiu^y, 1) 



(63) 
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and from the positive x region 



|^(0,2/) = Y2ik' h C h u h (y,u) 

oo 

= ik 'h c h u hh u h{vA)- 

h,h=o 

So then the continuity of dip/dx at x = can be written as 

( kA . . . ) - ( k B kxBx •••) = ( k' C Q Kd . . 



(^00 ^01 • • • 

v : ''•/ 



If we define the matrices 



k :-- 



k 
k x 



and 



( K o 



k' 



k[ 



\ 



J 



(64) 



(65) 



(66) 



(67) 



then we can write ( 1651) succinctly as AA; — B k = Ck'U. 

To solve the scattering problem (getting B and C in terms of A), we must solve ( 162]) and ( 165]) together. 
Solving for C, we get 

C = 2A(k)(Uk + k'Uy l . (68) 



Solving for B we get 



B = A(k- U- l k'U)(k + U^k'U)- 1 



(69) 



This gives us the S-matrix for this potential. 

We can analyze these solutions by computing the corresponding transmission and reflection coefficients. 
These are usually defined for one- dimensional scattering problems in terms of the probability flux in 
the x direction, which we can also do here. Recall that the probability current vector j is defined by 
ji '■= -y (ip*diip — ipdiip*) and its conservation is ensured by the Schrodinger equation. We can then use 
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FIG. 4. Plot of the partial reflection R4 for ui = 10.5 and matrices of size 10 x 10, 20 x 20, 40 x 40, 60 x 60, and 
80 x 80, shown in black, blue, purple, red, and gray respectively. 



Stoke's theorem on a region C bounded by some lines at x — ±c to obtain 

p poo /*oo 

= / V • j(x,y) dx dy = / j(-c, y) ■ (-x) dy + j(c, y) ■ (x) dy 

Jc J —00 J~ 00 

= -k\Af + ^Re(^)l^| 2 + ^Re(A;; i )|Q 1 | 2 . 



(70) 
(71) 



Identifying the total transmission coefficient 



T 



and the total reflection coefficient 



R 



E h Mk' h )\c h \ 

k\A\ 2 



k\A\ 2 



(72) 



(73) 



we have T + R — 1. We can also define the partial transmission and reflection coefficients for a given level 
I as Ti = Re(A; / , )|a| 2 A|A| 2 and Ri = Re(^)|^| 2 /A;|A| 2 . 

We can calculate these coefficients for a given incoming amplitude using ( )68|) and ( 169|) . but we run 
into the problem of inverting certain infinite matrices. We do this numerically by truncating the matrices 
to some finite size. We also ensure that the truncated version of the matrix U in ( 1531) is unitary by 
orthogonalizing (and normalizing) its columns via the Gram-Schmidt algorithm. This allows us to avoid 
numerical issues encountered because simply truncating U to finite size results in a matrix with very small 
eigenvalues. As the orthogonalized truncation is taken to large matrices, our results should converge to 
their true values. We have checked this numerically as shown in Figure HJ There we show the curve for 
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FIG. 5. Plot of the relative discrepancy (R^ with N x N truncation) /(JS4 with 80 x 80 truncation) for cj = 10.5 
and iV = 10, 20, 40, 60 and 80, shown in black, blue, purple, red, and gray respectively. 

the second partial reflection coefficient for u = 10.5 for matrices of increasing size. We see that the curve 
converges rather rapidly. To separate the curves better, we also show the relative discrepancy from the 
value calculated with 80 x 80 matrices in Figure [5) The lesson from numerous such numerical investigations 
is that, expectedly, one needs matrix sizes much larger than the highest channel one wants to accurately 
consider. 

In Figure |6] we display a logarithmic plot (base 10) of the first four partial transmission coefficients for 
the case A — 1, u — 2.5, for k values ranging from 1.0 to 10.0 and calculated using 50 x 50 matrices. On 
the horizontal axis we show the values of momenta k^ = uo(2l + 1) — 1 at which it becomes kinematically 
possible to transmit into channel I. As expected, this is where we see the corresponding partial transmission 
coefficient become non-zero. 

Also shown are the asymptotic values expected for these coefficients in the limit of large k. These come 
from the high k limit of ( |68|) . which is C = AU~ l . In this limit there is no reflection and it corresponds 
to the classical limit for x. This leads to T\ — > Uj^, as shown. This is the result one would expect from 
just performing a Bogoliubov transformation between the mode functions for a time dependent abrupt 
change in a potential. We expect these asymptotic values to decrease exponentially with level number (See 
Appendix 1). The proof that this is the correct result is that asymptotically we have that ki ~ k[ ~ ko for 
all modes, up to corrections of order 1/k. Thus Uk'U~ l ~ k and B n — > 0. 

Notice that at finite k values the conservation of energy plays an important role in forbidding transitions 
to high excitation. For these transitions the naive classical result of the Bogoliubov transformation gives 





FIG. 7. Logarithmic Plot of the first four partial reflection coefficients for the case u) — 2.5. Rq ^ R2 R4. Rq are 
shown in black, blue, purple, and red, respectively. 
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the wrong answer and the probabilities converge at different values of k for each mode. Our plots show 
how that works. The other thing to notice is that reflection coefficients are generally small. However, 
these would not be present at all in classical physics. 

The logarithmic plot of the partial reflection coefficients contains more structure, as shown in Figure 
[3 As expected the coefficients all approach zero in the limit of high k, corresponding to negative infinity 
on the plot. We have restricted the range of k to (1.4, 10.0) to increase the contrast. The vertical dashed 
lines correspond to the onset of transmission at momenta k\ as above. The vertical dotted lines indicate 
when a new reflection channel opens up, k* Rl = \f%. 




We would like to remind the reader that having reflection in this setup is a purely quantum effect and can 
not be obtained from a semiclassical expansion around a classical trajectory for x. Because we have solved 
the problem exactly, as given by equations (1681) and (1691) . our system has taken into account back-reaction 
to the quantum y motion exactly. We also can see that the system converges to the usual treatment of 
a time dependent harmonic oscillator in the classical limit for x and we can quantify the convergence 
precisely. The conservation of energy also tells us that each mode k[ will propagate at different speeds, so 
that the different components of the wave function will separate. This effect can be measured and further 
interactions with other degrees of freedom would probably decohere these different branches of the wave 
function. 

In a similar problem in field theory, the mass m x is related to the volume of space. The large mass 
limit is equivalent to the large volume limit. One can also state that the large mass limit is equivalent 
to a weak coupling limit by a rescaling of fields associated with x, but this does not mean much in this 
abrupt scenario. Notice that one gets finite particle production per unit volume and the general matching 
involves all the modes of the heavy fields. Taking the large volume limit is tricky because the x degree of 
freedom becomes part of a field. The right way to solve the problem is to coarse grain it into cells of fixed 
volume, so that gradients between adjacent cells are subleading in energy and so that the wavelengths 
of the particles produced are typically much shorter than the size of the cells. In such a way one can in 
principle perform the above calculation cell by cell and have the possibility of reflection or not on different 
boxes generating fluctuations in x. Notice that different positions will evolve afterwards at different rates 
(if one imposes energy conservation cell by cell), so the problem becomes rather involved even in the 
presence of a simplified abrupt potential as we have considered above. Such a calculation is beyond the 
scope of the present paper. 
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V. CONCLUSION 

From the simple system described by (pQ), for one modulus x and one harmonic oscillator (fast) degree 
of freedom whose frequency, Q(x), depends on x, we have gotten insight into the near-adiabatic and 
non-adiabatic dynamics of quantum systems. Following the method of Born and Oppenheimer we have 
been able to clearly see the transition from adiabatic to non-adiabatic behavior, including back-reaction, 
and to understand when the oscillator degrees of freedom cannot be integrated out. We implemented 
an improved adiabatic approximation that gave us a velocity dependent potential for the moduli motion. 
This accounted for some back-reaction effects. We saw how this back-reaction depended on the velocity 
of the modulus x and the coupling between x and the oscillator, as well as how it can be interpreted as 
modifying the metric on the moduli space. We also argued that for finite time evolutions one obtains 
an occupation number for the oscillator degrees of freedom that is a power series in the non-adiabaticity 
parameter. These effects are responsible for the velocity dependent corrections to the Hamiltonian when 
the oscillator degree of freedom is integrated out. However, for infinite time behavior the net particle 
production can be non-perturbatively small. 

Despite the simplicity of the model, which we have analyzed in a non-relativistic classical and quantum 
mechanical setting, it is general enough to be relevant to both relativistic quantum field theory and string 
theory as we have briefly discussed and as evidenced by its frequent appearance (for various choices of Q) 
in the literature. We plan to apply these results to these setups in the future. In particular, we found 
that the velocity corrections have a definite sign, the same for fermions and bosons. However, we did not 
carry out a complete analysis for gauged systems, which may be behave differently in important ways 
(e.g. because of the presence of ghosts). This is especially important to address because of the non- 
renormalization theorems that constrain the dynamics, via cancellations, when supersymmetry is present 



23]. It is important to understand exactly how these cancellations happen in supersymmetric gauged 
systems, as in those systems associated with D-branes. 

We have also solved a fully non-adiabatic problem exactly and investigated the solution numerically. 
Our model had an abrupt change in Q that gave rise to quantum mechanical transmission and reflection, as 
in one-dimensional scattering problems, which are displayed in the final four figures. The kinematics gave 
rise to a separation of different modes, leading to possible decoherence of the wavefunction amplitudes as 
time goes on and the system interacts with other degrees of freedom. It is easy from our solution to treat 
the interesting case of an Q with multiple abrupt changes and we plan to investigate this in future work. 

Many of the phenomena we have encountered are inherently quantum-mechanical and may be important 
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to consider in applications, such as D-brane interactions or inflationary cosmology, where previous analyses 
have been semiclassical in nature. It may be crucial in making a problem tractable that these effects are 
suppressed, so that one may get by with semiclassical physics. We think it is thus important to examine 
very carefully when and to what extent this is the case. We found the above model, which admitted exact 
quantum solutions in certain simple situations, an illuminating test case for such an investigation. 
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VI. APPENDIX 1 



The Hamiltonian for the y variable in the two regions separated by the y-axSs differ only in their value 
of Q. Let us call the region where x < region 1 and where x > region 2 and consider a Hamiltonian 
with two different values for Q, oj\ and u 2 respectively. Since we are interested just in ensuring continuity 
of the wavefunction at x = we can ignore the x dependence of the states here. The problem then is 
to find the transformation that will take us from the basis of energy eigenstates, which we can write as 
{|n)i}, in one region to the basis {1^)2} m the other. That is, we want to calculate the matrix in ( !53|) : 

U mn = (n\im) 2 = ^=(0| 1 (a 1 )™(4) m |0) 2 . (74) 
\/n\m\ 

This can be accomplished by relating the creation and annihilation operators as well as the vacuum states 
on the two sides. For i — 1, 2 we have 

~~ y + ^) (75) 
4 = J^(v-^)- P6) 





Then we can calculate 




Ul + u 2 oo\ - u 2 



V^i^s \ Ul - u) 2 Ui + UJ 2 




(77) 
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a Bogoliubov transformation |4|, |21[ or quasiparticle transformation |20{|. To relate the vacuum states 



1 0), we note that they are squeezed states with respect to the opposite operator algebra. To find an 
explicit formula for |0)i we recall the representation of states by holomorphic functions [20j and write 
1 0) 1 = F(a\)\0) 2 for F holomorphic. Then since formally a 2 = d/da) 2 when acting on the vacuum state, 
ai|0)i = yields the equation for F 



d 

(ui + u 2 ) — + (wi -u 2 )z) F(z) 



0. 



(7* 



which has solution F{z) = C exp ^- 



1 u>i~u) 2 ^2 

2 u>i+ui2 



Normalizing, we find 



V Ui + tu 2 J 



1/2 



exp - 



1 Uli — UJ 2 



2ui+u 2 



|0) ; 



(79) 



Thus 



1/2 



Un 



1 / 2y/U^Up \ 

y/n\m\ \ ^1 + ^2 / 

1 U)\ 



x ( 



exp 



2wi + w 2 



U 2 t 2 

a2> 



[ui + u 2 )a 2 + (ui - 0J 2 )a\ 







10) 



This expresssion can in principle be evaluated using the commutation relations for a 2 and a 2 , though we 
found it easier to work directly with the Hermite polynomial generating function to derive ( 15^1) . But (IHUl) 
is useful in determining the large m or large n asymptotics of the off-diagonal elements and hence of the 
partial transmission coefficients considered in section IIVI 

First we note a certain symmetry property of the matrix U mn , which is easiest to see from the inte- 
gral expression (|53l) . Since the functions H n (x) exp(— x 2 /2) are eigenfunctions of the (unitary) Fourier 
transform, with eigenvalue (— i) n , Fourier transforming the functions in ( l53p yields 

1 



Un 



—I 



im+n 



2tt 




u n (p)u m (p')e ix{p+p '^ dx dp dp' 



—i 



•.rn+n 



(81) 



This says U mn is symmetric up to a phase. 

Since here we are just interested in the absolute value of the off-diagonal elements, by the above 
symmetry property we lose no generality if we fix n and consider large m. The leading contribution comes 
from the term in the expansion of the exponential in f lHUj) with the appropriate power to annihilate the aj> 
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operators. We take m 3> n and focus on the m-dependent factors, using the Stirling approximation on all 
factorials. We find that the m dependence is asymptotically 

which decreases exponentially in magnitude with m. 
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